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ABSTRACT: 

A new type of monolithic optical waveguide filter comprises a chain of optical couplers (14,15,16) 
of different effective lengths linked by differential delays (17,18) of different lengths. It can be 
regarded as a series of Mach - Zehnder interferometers. The transfer function of the chain of 
couplers and delays is the sum of contributions from all possible optical paths, each contribution 
forming a term in a Fourier series whose sum forms the optical output. A desired frequency 
response is obtained by optimizing the lengths of the couplers and the delay paths so that the 
Fourier series best approximates the desired response. The filter is advantageously optimized so 
that it is insensitive to uncontrolled fabrication errors and is short in length. The wavelength 
dependence of practical waveguide properties is advantageously incorporated in the 
optimization. Consequently, the filter is highly manufacturable by mass production. Such filters 
have been shown to meet the requirements for separating the 1 .3 and 1 .55 mu m 
telecommunications channels and for flattening the gain of Er amplifiers. 
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(54) Monlithic optical waveguide filters based on Fourier expansion 



(57) A new type of monolithic optical waveguide fil- 
ter comprises a chain of optical couplers (14,15,16) of 
different effective lengths linked by differential delays 
(17,18) of different lengths. It can be regarded as a 
series of Mach - Zehnder interferometers. The transfer 
function of the chain of couplers and delays is the sum 
of contributions from all possible optical paths, each 
contribution forming a term in a Fourier series whose 
sum forms the optical output. A desired frequency 
response is obtained by optimizing the lengths of the 
couplers and the delay paths so that the Fourier series 



best approximates the desired response. The filter is 
advantageously optimized so that it is insensitive to 
uncontrolled fabrication errors and is short in length. 
The wavelength dependence of practical waveguide 
properties is advantageously incorporated in the optimi- 
zation. Consequently, the filter is highly manufacturabJe 
by mass production. Such filters have been shown to 
meet the requirements for separating the 1.3 and 
1.55^im telecommunications channels and for flattening 
the gain of Er amplifiers. 
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Description 

Field of the Invention 

This invention relates to monolithic optical 
waveguide titters and, in particular, to a new type of 
monolithic fitter providing plural optical paths wherein 
each optical path corresponds to a harmonic compo- 
nent of a Fourier series comprising the filter transmis- 
sion function. 

Background of the Invention 

Optical filters are important devices in optical ffoer 
communications systems. Monolithic optical waveguide 
filters are particularly promising because they can per- 
form complex circuit functionalities and because they 
can be made by mass production integrated circuit tech- 
niques. 

The requirements of optical filters vary with applica- 
tions. Many applications require a rectangular wave- 
length response in order to maintain a low-loss and 
wavelength-independent transmission in a passband 
and a high-level rejection to all wavelengths in a stop- 
band. For example, anticipated telecommunications 
applications seek a 1 .3/1 .55^m WDM filter a flat and 
low-loss passband at 1.280-1.335 (am and a -50 dB 
stopband at 1.525-1.575 fim. Another desired applica- 
tion is a gain equalization filter to flatten the gain of an 
Er<Joped fiber amplifier chain. This requires an equali- 
zation filter with an amplitude response which is essen- 
tially the inverse of the amplifier gain. 

Various devices have been proposed to fill these 
new, demanding requirements but none are fully satis- 
factory. Multilayer thin-film filters can be used to con- 
struct optical filters in bulk optics, but they are 
undesirable because they cannot be readily integrated 
and because of difficulties in coupling light to and from 
fibers. Mach-Zehnder (MZ) interferometers have been 
widely employed, but they have a sinusoidal response, 
giving rise to strongly wavelength-dependent transmis- 
sion and a narrow rejection band. Other designs have 
encountered a variety of practical problems. Accord- 
ingly, there is a need for a new type of monolithic optical 
waveguide filter. 

Summary of the Invention 

A new type of monolithic optical waveguide filter 
comprises a chain of optical couplers of different effec- 
tive lengths linked by different differential delays. The 
transfer function of the chain of couplers and delays is 
the sum of contributions from all possible optical paths, 
each contribution forming a term in a Fourier series 
whose sum forms the optical output. A desired fre- 
quency response is obtained by optimizing the lengths 
of the couplers and the delay paths so that the Fourier 
series best approximates the desired response. The fil- 
ter is advantageously optimized so that it is insensitive 



to uncontrolled fabrication errors and is short in length. 
The wavelength dependence of practical waveguide 
properties is advantageously incorporated in the optimi- 
zation. Consequently, the filter is highly manufacturable 
5 by mass production. Such filters have been shown to 
meet the requirements for separating the 1 .3 and 1.55 
telecommunications channels and for flattening the 
gain of Er amplifiers. 

10 Brief Description of the Drawings 

In the drawings: 

FIGs. 1 a, 1 b and 1 c are schematic design layouts of 
15 1 .3/1 .55 p.m wavelength division multiplexing filters 
(WDM filters); 

FIGs. 2a and 2b are diagrams useful in illustrating 
the design principle that the filter output is the sum 
of all optical paths; 
20 FIGs. 3a, 3b and 3c are schematic examples of 
four-coupler chain filters of consecutive odd Fourier 
harmonics; 

FIGs. 4a, 4b, and 4c illustrate the approximation of 
a rectangular fitter response by a Fourier series; 

25 FIGs. 5a, 5b, 5c, 5d and 5f show the effect of vari- 
ous steps to optimize a five coupler chain filter; 
FIG. 6 illustrates the basic functions of a coupler- 
delay chain as a four port optical filter; 
FIGs. 7a, 7b and 7c are graphical illustrations of the 

30 transmission spectra for embodiments of WDM fil- 
ters of the designs shown in FIGs. 1a, 1b and 1c, 
respectively; 

FIG. 8 shows the amplitude response of a gain 
equalization filter to flatten the overall gain of Er- 
35 doped amplifiers; and 

FIG. 9 is a schematic top view of a simple fitter. 

Detailed Description 

40 This description is divided into five parts. In part I, 
we describe the basic elements of a simple filter. In part 

II we describe the physical fabrication of the filter. Part 

III is directed to the design of the configuration of 
waveguides to obtain a desired fitter response. Part IV 

45 discusses practical considerations which assist in the 
fabrication of practical filters; and Part V discusses pre- 
ferred uses of the filters in optical fiber communications 
systems. 

so I The Basic Elements Of A Simple Filter 

Referring to the drawings, FIG. 9 is a schematic top 
view of a simple form of a monolithic optical waveguide 
filter 10 comprising a pair of optical waveguides 1 1 and 
55 12 on a substrate 13 configured to form a plurality N of 
at least three optical couplers 14, 15, and 16 alternately 
connected by a plurality of N-1 delay paths 17 and 18. 
Each coupler is comprised of a region of close adja- 
cency of the two waveguides where the exponential tail 
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of light transmitted on each of waveguides 11 and 12 
interacts with the other, coupling light from one 
waveguide to the other. The amount of power coupled 
from one waveguide to the other is characterized by the 
effective length of the coupler. The effective lengths of 
the couplers preferably differ from each other by more 
than 5%. The effective length of any coupler is within 
5% of at most one other coupler. 

Each delay path comprises a pair of waveguide 
segments between two couplers, for example segments 
17A and 17B between couplers 14 and 15. The seg- 
ments are configured to provide unequal optical path 
lengths between the two couplers, thereby providing a 
differential delay. For example in FIG. 9 upper segment 
17A is longer than lower segment 17B, providing a dif- 
ferential delay which can be denoted positive. Differen- 
tial delays associated with longer lower segments can 
be denoted negative. Delay path 18 provides a negative 
delay because lower segment 18B is longer than upper 
segment 18A. In the preferred form of filter 10, at least 
one differential delay differs from at least one other by 
1 0% or more and at least one differential delay is oppo- 
site in sign from at least one other. 

In operation, an optical input signal is presented at 
an input coupler, e.g. along waveguide 1 1 to coupler 14, 
and a filtered output is presented at an output coupler, 
e.g. along waveguide 1 2 at coupler 1 6. The sequence of 
couplers and delays provide light at the input with a plu- 
rality of paths to the output. In general there will be 2 N " 1 
paths where N is the number of couplers. For example, 
the FIG. 9 device presents the following four paths: 

1} segment 17A - segment 18A 

2) segment 17A - segment 18B 

3) segment 17B - segment 18A 

4) segment 17B - segment 18B 

Each of the optical paths of the filter provide light 
corresponding to a harmonic component in a Fourier 
series whose summation constitutes the transmission 
function of the filter. By proper choice of parameters one 
can closely approximate a desired transmission func- 
tion. By choice of N and the set of differential delays, 
one can design a filter presenting a summation of odd 
Fourier components particularly useful for fabricating a 
fitter with a rectangular response such as a 1 .3/1 .55 jim 
WDM filter. With a different N or a different set of differ- 
ential delays, one can also design a filter presenting a 
summation of both even and odd Fourier components. 
A preferred odd harmonic filter can be made by provid- 
ing differential delays with a normalized ratio 
±1£2£2/- • "h2. in any order, and a preferred all har- 
monic filter can be made by providing differential delays 
with a normalized ratio of ±1#1/+2/±2/- • */t2 in any 
order. Moreover, in the above ratios, any but not all of 
the differential delays of ±2 can be replaced by ± 4 and 
any but not all of the differential delays of ±4 can in turn 
be replaced by ±8. In addition, for broadband filters, the 
wavelength dependence of practical waveguide proper- 



ties is advantageously taken into account, which can 
alter the above proportional differential delays by up to 
±25%. Thus in one preferred embodiment the delay 
paths provide, in any order, within ±25%, one differential 

5 delay of proportion ±1 and one or more differential 
delays of proportion ±2 or ±4 or ±8. In another preferred 
embodiment the delay paths provide, in any order, 
within ±25%, one differential delay of proportion ±1 , one 
more differential delay of proportion ±1, and one or 

10 more cfifferential delays of proportion ±2 or ±4 or ±8. 
Advantageously the filter can be combined with other fil- 
ters, as by connecting the other filters to the waveguide 
outputs, thereby producing filter networks. 

15 II. Physical Fabrication 

The FIG. 9 structure is advantageously fabricated 
using planar optical waveguide technologies. Doped sil- 
ica optical waveguides are preferred because they have 

20 low loss, low birefringence, are stable, and can couple 
to standard fibers well. However, the invention is equally 
applicable to other integrated optical waveguides 
including lll-V semiconductor optical waveguides and 
optical waveguides diffused in lithium niobate. A 

25 description of the above waveguide technologies can be 
found in Ft. G. Hunsperger, "Integrated Optics: Theory 
and Technology", 3rd ed. (Springer-Verlag, Berlin, Hei- 
delberg, New York 1991). 

With doped silica plannar waveguides, the FIG. 9 

30 structure can be fabricated much as described in C. H. 
Henry et al., "Glass Waveguides On Silicon For Hybrid 
Optical Packaging," Lightwave Technol. . vol. 7, pp. 
1530-39 (1989). In essence a base layer of silica glass 
(SiO^ is grown on a silicon or quartz substrate. A thin 

35 core layer of doped silica glass is then deposited on the 
base layer. The core layer can be configured to a 
desired waveguide structure, such as that shown in FIG. 
9, using standard photolithographic techniques. Subse- 
quently another layer of silica glass is deposited to act 

40 as a top cladding. The waveguide cores have a higher 
refractive index than the base and top cladding layers, 
thereby guiding the lightwave much as a fiber. In con- 
trast to a fiber, however, planar waveguide technologies 
are capable of more precise control of the waveguide 

45 structures and of integrating many devices into complex 
circuits. 

In the specific examples discussed below, the fol- 
lowing procedure was used to fabricate the filter. First a 
~15*un thick base layer of undoped Si0 2 (HiPOX) is 

so formed by oxidation of Si under high pressure steam. A 
core layer of 5 \im thick 7% P-doped Si0 2 (p-glass) is 
then deposited using low-pressure chemical vapor dep- 
osition (LPCVD). The core layer is annealed in steam at 
1000° C is patterned tor single mode waveguides of 5 

55 jim width with reactive ion etching (RIE). The center-to- 
center separation of waveguides in the coupler regions 
is 9.25-9.5fim, and the minimum waveguide bend radius 
is 15 cm. A flow layer of about 1.5 jwn thick B- and P- 
doped Si0 2 (BPTEOS) is subsequently deposited with 
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LPCVD and annealed at 900° C to fill the narrow gaps 
in between waveguide cores in the coupler regions. 
Finally, two layers of BPTEOS of 7.5 \im thickness each 
are deposited and annealed as the top cladding. The 
top cladding has almost the same refractive index (1 .45 
at X=1.4fxm) as the base, and the core has a refractive 
index of about 0.63% higher than that of the base and 
top cladding. 

ill Confi guration Design And Optimization 

The basic structure of our filters consists of a chain 
of N arbitrary couplers and N-1 differentia! delays, 
where N>2. The transmission spectrum of such a chain 
(referred to as an N-coupler chain) is the sum of contri- 
butions from 2 N1 optical paths, each of which forms a 
term in a Fourier series. The length of the couplers and 
delay paths can be optimized so that this Fourier series 
best approximates a desired frequency response. As an 
initial example we describe design of the aforemen- 
tioned 1.3/1.55 |xm WDM filter with a rectangular 
response. FIGs. 1a, 1b and 1c are layout examples of 
1 .3/1 .5 jim WDM filters. FIG. 1 a shows a basic five cou- 
pler chain. Arrows indicate input and output ports. FIG. 
1b shows a doubly filtered five-coupler configuration. 
FIG. 1c show a combination of double and triple filtering 
of four-coupler and five-coupler chains. When fabri- 
cated with the aforementioned doped silica waveguide 
technology, the total lengths of the 1a, 1b and 1c filters 
are 43, 75 and 75 fim, respectively, and the heights are 
0.3, 0.6 and 0.6 mm, respectively. The vertical scale has 
been expanded 20 times for clarity. 

The Principle nf Sum Over All Possibl e Optical Paths 

The electric fields at the two output ports of an ideal 
coupler of two identical waveguides are related to those 
at the input ports by a transfer matrix 



rcos <(> i sin (H 
u sin $ cos <|>J 



(1) 



where i=V-1 . ♦ = rcl72L,r is the geometric length of the 
coupler, and L is the coupling length which is a measure 
of the strength of coupling between the two waveguides 
(not to be confused with the geometric length of the 
coupler). We refer to I'/L as the effective length of the 
coupler. 

Similarly, the transfer matrix characterizing the dif- 
ferential delay between two identical waveguides is 



0 e 



(2) 



where e=7csnA,=Tisnv/c , s is the difference in the 
lengths of the two waveguides, n is the effective refrac- 
tive index of the waveguides, and X, v, and c are the 
optical wavelength, frequency, and velocity in free 



space. Note that a common phase factor has been 
ignored in Eq. 2 because it is non-essential to the filter 
response. 

The transfer matrix of our waveguide filter, consist- 
ing of a chain of N couplers and N-1 differential delays 
characterized by <h.te." * •.♦n. and e i- 0 2- * " ". e N-i- 
respectively, is given by 



T = T T 

4>o 4> N (Jn-i 



■T H T A T 0 T. 

B 2 92 U 1 4*1 



(3) 



From the above matrix product we see that the transfer 
function from any input port to any output port consists 
of a sum of the form 

v K±o,±e 2 ' • -±0^) 
* t +0 = Lf(ti.*2.---^N)e < 4 > 

Note that 6 is proportional to v which is proportional to 
MX if we ignore the wavelength dependence of n. 

Physically, equation 4 can be interpreted as the fol- 
20 lowing principle that the response is the sum over all 
possible optical paths. The transfer function from any 
input port to any output port of a chain of N couplers and 
N-1 differential delays consists of the unweighted sum 
of contributions of 2 N1 distinct optical paths. Each of 
25 such contributions is a product of 2 N-1 factors: travers- 
ing a coupler gives cos<t> without crossing and isin<|> with 
crossing; traversing the longer arm of a differential delay 
gives e* 6 and the shorter arm gives e"'° 

This principle is illustrated in FIGs.2 (a) and 2 (b). 
30 FIG. 2a diagrams the contributions from the paths in a 
coupler and a differential delay. FIG. 2b shows the cross 
state of a three coupler chain. Where we have used the 
abbreviations c«cos<t> and s=sin<|> . The letter t desig- 
nates the transfer function. For N=3 there are four dis- 
ss tinct optical paths from any input port to any output port, 
and the transfer function shown in FIG. 2(b) is a sum of 
four terms. For non-ideal couplers, cos<|> and isirkj) 
should be corrected accordingly, but the above principle 
is still valid. 

40 

Construction of a Fourier Series 

With arbitrary choices of the 6's, the sum in Eq. (4) 
is normally not a Fourier series because the terms in the 
45 sum do not normally represent harmonics. However, 
Eq. (4) becomes a truncated Fourier series if the ratio 
among the lengths of the differential delays satisfies 
certain conditions. In the following, we define 



so 



ejM-YjeoMJ-IA- • - ,N-1, 



(5) 



where % contains the common wavelength-dependent 
part of 8-1,62, • • • ,e N _-| and has the same period as the 
fundamental harmonic in the Fourier series. The unnor- 
55 malized ratio y^/ • • • /yn-1 will be referred to as the 9- 
ratio. 

For WDMs with a rectangular response, because of 
the odd symmetry of the required transfer function, we 
need a Fourier series with only odd harmonics (see 
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FIG. 4). We first consider the following two extreme con- 
ditions which give consecutive odd harmonics: 

Condition A : K N is even and the 0-ratio is 
±1y±1/- * *l± 1 with any sign combinations, i ¥) is a Fou- 
rier series with harmonics 
±0 O ±36 0l • • -,±^-1)60. 

Condition B : If the 0-ratio is ± 1/ ±2/ ±4/- • -/±2 N - 2 , 
in any order and with any sign combinations, t^o is a 
Fourier series with harmonics 
±e 0 ,±3eo • • • f ±(2 N - 1 -1)0 o . 

For a given even N, condition B gives the maximum 
number of consecutive odd harmonics while condition A 
gives the minimum. There are many other conditions in 
between, such as: 

Condition C : If the 0-ratio is ±1/ ±2/ ±2/ / ±2, 

in any order and with any sign combinations, t^) is a 
Fourier series with harmonics 
±e 0i ±36 0( • • -,±(2N-3)0 O . 

FIGs. 3a, 3b and 3c illustrate the construction of the 
Fourier harmonics using N=4 with examples of 0-ratios 
of 1/1/1 , 1/2/4, and 2/1/-2, each corresponding to a spe- 
cial case of conditions A, B, and C, respectively. Note 
that a negative sign in the 0-ratio corresponds to an 
interchange of the longer and shorter delay arm. For 
clarity, we have used 0 in the drawing for 0 O . A negative 
<|> cannot be physically realized. However, is equiva- 
lent to in that the transfer functions of the two cases 
differ only by a constant phase. 

While consecutive odd harmonics are needed for 
WDMs, other filter functions may need even harmonics 
as well. A Fourier series with consecutive even harmon- 
ics is constructed when one of the sections correspond- 
ing to the ±1 's in conditions A, B, and C is taken out of 
the N-coupler chain, resulting in an (N-l)-coupler chain. 
A Fourier series with all consecutive harmonics can also 
be constructed by halving all the delays in conditions A, 
B, and C and adding another section with the minimum 
unit delay. The 0-ratios for conditions A, B, and C 
become ±^±W* ■ */± 1 ^ (and N is odd), ±V4/±V4/±1/ 
±2/ — /±2^" 4 ,and± 1 >4/± 1 >4/±1/ • • WH, respectively. 
The corresponding normalized ratios of the differential 
delays are ±1/±1/* • */±1 ,±1/±1 1+21 ±4h • -/±2 N - 3 , 
and+1/±1/±2/* • */±2, respectively. 

Optimization for a Given Filter Response 

We have shown that the transfer function of a chain 
couplers and delays can form a truncated Fourier 
series. FIG. 4 further demonstrates the basis of approx- 
imating the frequency response of a desired filter by 
such a truncated Fourier series (for the purpose of illus- 
tration, we have ignored the phase). The MZ interferom- 
eter (N = 2) only has the fundamental harmonic and its 
frequency response sinusoidal as shown by the solid 
curve in FIG. 4(a). In contrast, the desired rectangular 
response of a filter is shown by the solid curve in FIG. 
4(b). For N > 2 under conditions A, B, or C, the transfer 
function also contains higher order harmonics, shown 
by the broken curves in FIG. 4(a). The sum of the funda- 



mental and higher order harmonics, as shown by the 
broken curve in FIG. 4(b), can approximate the rectan- 
gular response if their amplitudes and phases are cho- 
sen correctly. In the following, in order to give a clear 

5 physical picture, we first ignore the wavelength depend- 
ence of L (and hence the <t>'s are constants) and n and 
discuss the more general case later. 

First, the common factor 0 O of the differential delays 
can be determined solely by the positions of the pass- 

to band and the stopband, similar to the case of a simple 
MZ. We want the fundamental harmonic to have a 
phase of 0 or n at the center of the stopband and k/2 at 
the center of the passband. Therefore 

75 0^/71 -sn/X s =m, and 0 Op /n-snA p =m±V£, (6) 

where m is an integer, which we refer to as the order of 
the filter, and the subscripts p and s refer to the center 
of passband and stopband. Since the output ports can 

20 be interchanged, there is another configuration corre- 
sponding to the interchange of s and p in Eq. 6. For our 
1 .3/1 .55 fim WDM filter, the best solutions are (m = 3, Ap 
= 1.322 jim, ^=1.542^) and (m = 3, Ap = 1.566fxm, A*. 
= 1 305^m). Here we have explicitly used the fact that 

25 the transfer function expressed by the Fourier series is 
periodic in frequency. 

We will now find the values of the <t>'s under a given 
condition that renders the Fourier series with consecu- 
tive odd harmonics. Note that, except under condition A, 

30 the number of harmonics exceeds the number of cou- 
plers N. Therefore, the number of Fourier coefficients to 
be determined is generally larger than the number of 
free variables available. Moreover, for broadband filters, 
the phase response is not important, and we only 

35 require It^ol 2 to have the desired response which is rec- 
tangular in our 1.3/1.55 WDM. We therefore minimize 
an error function to solve for the <t>*s under a given 0-ratio 
as in the following: 

« E^-/dv * w(v)(|t^|-t desirGd (v)|) 2 =min (7) 

where w(v) is a positive weighting function, and the inte- 
gration is done in the passband and stopband of inter- 
est. Since for our 1 .3/1 .55fim WDM we do not constrain 

45 the transition between the passband and stopband, we 
set w(v)=0 in those regions (as in FIG. 4(c)). Eq.(7) rep- 
resents a nonlinear minimization problem which can be 
solved numerically by iterative methods such as the 
simplex or conjugate gradient methods. Although in our 

so design of the 1.3/1.55 \im WDM we have ignored the 
phase response, if a particular phase response is 
desired, it can also be put into the above equation and 
optimized together with the amplitude. 

Note that E^ has many local and equivalent 

55 minima which correspond to different configurations 
having similar or the same amplitude response. For 
example, if one of the couplers in the middle of an N- 
coupler chain is a full coupler (ty={n+ ] £)K where n is an 
integer) or a null coupler (<j>=n7x ), this coupler and the 
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two differential delays around it degenerate into the 
equivalence of a single differential delay, and we effec- 
tively have an (N-l)-coupler chain. Thus we can repro- 
duce a minimum in E^, of an (N-1 )-coupler chain with an 
N-coupler chain, which is a local minimum. Moreover, rf 
the length of a coupler is changed in such a way that $ 
is replaced by <t>+2nn or (2n-1)7nj>, then \t^\ is 
unchanged, and we have an equivalent minimum. Fur- 
thermore, a different sequence of 9's or a different set of 
± signs in front of the e's gives different equivalent and 
local minima. The different combinations of these varia- 
tions give rise to large number of configurations. Inter- 
changing the two output ports brings in another set of 
distinct configurations (which correspond to $^n/4 
and <t> 2 ^ 37t/4 for the conventional MZ). Therefore, 
physical insight is often needed to get good starting 
parameters for the iterations to arrive to the best solu- 
tion. We usually start with the original MZ (N = 2) and 
successively increase N, taking the results of N - 1 as 
the starting parameter of N. 

We still have to find the best choice of the e-ratio 
that gives the closest approximation of the required filter 
response for a given N. In the extreme of condition A, all 
the Fourier coefficients can be satisfied independently, 
but the number of harmonics in the series is small. In 
the other extreme, under condition B, although the 
number of harmonics is large, they are greatly con- 
strained by the number of free variables. Therefore, 
there is an optimum condition in between A and B. Our 
numerical calculations show that condition C is in fact 
the optimum for filters with a rectangular response. 

FIGs. 5a through 5f illustrate the effect of various 
steps to optimize the configuration design. FIG. 5a 
shows the bar and cross state power transmission of a 
conventional Mach-Zehnder interferometer in a semi- 
log plot. 

The solid curves in FIG. 5(b) are the bar- and cross- 
state transmission spectrum of a five-coupler chain sat- 
isfying condition C with a 9-ratio of M21-2I-2 (see FIG. 
1(a) for the layout), optimized when assuming L and n 
are independent of X. Compared to the spectrum of the 
corresponding MZ shown in FIG. 5(a), the width and 
flatness of the passband, as well as the width and rejec- 
tion of the stopband, are greatly improved. The dashed 
curves in FIG. 5(b) show the corresponding transmis- 
sion of a six-coupler chain with a 9-ratio of 1/1/1/1/1 
(i.e., condition A). Compared to the solid curves of the 
five-coupler 1/2/-2/-2 chain, the passband becomes 
narrower and the rejection lower. In fact, the response of 
the this six-coupler chain is the same as a four-coupler 
chain with a 6-ratio of 1/2/2. Similarly, when we replace 
one or more of the 2's in the 1/2/-2/-2 chain by 4 or 8, we 
find that the filter response becomes worse as well as 
the device length becomes longer. We therefore con- 
clude that condition C gives the most efficient WDM 
configuration. 



Wavelength Dependence of the Coupling Length and 
the Effective Refractive Index 

So far we have assumed that the coupling length L 

5 and the effective refractive index n do not depend on 
wavelength, which is only valid for narrow-band WDMs. 
In reality, L inevitably decreases as the the wavelength 
is increased, because the optical field is more confined 
at shorter wavelength. With our planar waveguide fabri- 

10 cation process, L increases by about a factor of two as 
X changes from 1.55fim to 1.3u.m, implying that a 3dB 
coupler at 1.3^im becomes a full coupler at about 
1.55nm. Moreover, n also changes with X because of 
the change in confinement and the dispersion of the 

15 waveguide material, but the relative change is much 
smaller (about 0.5 % for our waveguide) than L FIG. 
5(c) shows the response of the same five-coupler chain 
as in FIG. 5(b) (which was designed for 
L = L k=i.42 H m = constant and n = n| Xm ^ A2VLm = con- 

20 stant), except now the wavelength dependent L and n 
are used in calculating the spectrum. It is apparent that 
for wide-band WDMs, such as our 1.3/1.55jim WDM, 
the X dependence of L and n is advantageously consid- 
ered in the design. 

25 With the X dependent L and n, each term in Eq. (4) 
is no longer completely periodic, and it becomes neces- 
sary to optimize the e's in Eq. 6 as well as the <|>'s to get 
the best filter response. This is illustrated by the solid 
curves in FIGs. 5(d) and (e). FIG. 5(d) shows the 

30 response of the five-coupler 1/2/-2/-2 chain when only 
the coupler lengths have been optimized with the X- 
dependent L and n, while FIG. 5(e) shows the same 
when both the coupler lengths and the delay lengths 
have been optimized. Only in the latter case have we 

35 recovered the rectangular response comparable to that 
for constant L and n. With our planar waveguide for the 
1.3/1.55|im WDM example, the e-ratio change is from 
1/2/-2/-2 to 1.187/1 .978M.849/-2.031. In the re-optimi- 
zation, we usually use the results for constant L and n 

40 as starting parameters. 

In summary, in our optical filter with a chain of arbi- 
trary couplers and differential delays, the basic building 
block has a quasi-periodic transfer function. The trans- 
fer function of the chain is the sum of contributions from 

45 all possible optical paths, each of which can form a term 
in a Fourier series. The task of designing a filter is to 
optimize the lengths of the couplers and the differential 
delays so that this Fourier series best approximates the 
desired filter response. Fourier expansion not only gives 

so a direct and intuitive description of the physical principle 
of the device, but also provides powerful and flexible 
design procedures. Our Fourier expansion approach 
enables us to find the most efficient e-ratio (which is not 
1/1/1/1/- • • used in lattice filters), to include the X- 

55 dependence of L and n which renders the e-ratio non- 
integral, and, as will be shown next, to arrive at fully opti- 
mized configurations most immune to dominant fabrica- 
tion errors which involve negative signs in the e-ratio. 
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IV. Practical Considerations 

Because of the dependence of L and n, many of 
the equivalent configurations discussed earlier become 
inequivalent, i.e., the degeneracies are removed. We 
thus have many filter configurations which give some- 
what different response. The number of such configura- 
tions is large. 

The transmission of two of these configurations for 
our 1.3/1.55fun WDM is shown in FIGs. 5(e) and (f). 
FIG. 5(f) shows a five-coupler chain the same as that in 
FiG. 5(e) except the 6-ratio is approximately 1/2/2/2. 
The filter response is also similar to that in FIG. 5(e) 
except for some fine details. However, a dramatic differ- 
ence develops between the two cases when the cou- 
pling length deviates from the nominal values, as shown 
by the dashed curves in FIG. 5(e) and (f) for which a 
10% overall increase in L has been used in calculating 
the spectra. The design in FIG. 5(e) is less sensitive to 
the change in L than the design in FIG. 5(f). For practi- 
cal fabrication, the coupling length is the parameter 
most susceptible to uncontrolled fabrication errors (e.g., 
errors in the geometry and refractive index of the core), 
and the design in FIG. 5(f) would have a low yield. 
Therefore, we always choose out of many configura- 
tions the few which are most stable against overall 
changes in the coupling length as well as with the best 
filter response and short couplers. 

We have designed filters with various response 
shapes and bandwidths, and for all cases the configura- 
tions most stable against overall changes in L have neg- 
ative signs in the 8-ratio. While the problem of stability is 
complicated in nature partly due to the wavelength- 
dependence of L, a simple physical interpretation is as 
follows. For a filter with rectangular response to be sta- 
ble against fabrication errors, the partial sum of the fun- 
damental and, successively, those of the low-order 
harmonics in the Fourier series of the transfer function 
should be stable near their zeros, ff half of the 9's have 
negative signs, the fundamental and low-order harmon- 
ics correspond to the optical paths crossing the smallest 
number of couplers, which can be regarded as the dom- 
inant optical paths. The zeros of the partial sums pro- 
duced by these dominant paths are stable against 
errors in L because they tend to depend only on the 
ratio of the lengths of the small number of crossed cou- 
plers. Filters with negative 6's usually also have short 
couplers. For the stable 1.3/1.55jim WDM example 
shown earlier, the e-ratio is 1.187/1.978/-1.849/-2.031. 
This stable design is important for manufacturability. 

Our design approach based on Fourier expansion 
also gives us clear guidelines to determine the number 
of stages of the chain. In principle, as the number of 
coupler stages is increased, the flatness of the pass- 
band and the rejection in the stopband will be constantly 
improved. In practice, however, several factors limit N 
from being too large. The first limitation is chip size. With 
our current fabrication process on a five-inch wafer, the 
largest N is about 14 for short delay arms (such as 



those used in our 1.3/1 ,55\im WDMs) without using 
waveguide U-turns (which would consume a large 
space on the wafer). This limit of N can be increased, 
however, if high-delta waveguides or U-turns are used. 

5 The second limitation is excess insertion loss, due to 
the bends at the ends of each coupler and intrinsic loss 
in the waveguides. A more subtle yet important consid- 
eration is fabrication accuracy and non-ideal effects of 
couplers and delays. As N is increased, higher order 

10 Fourier components, with smaller Fourier coefficients, 
come into play However, rf the fabrication error exceeds 
the accuracy required by the smallest Fourier coeffi- 
cient, increasing N no longer improves the performance 
of a real device. Similarly, because of non-ideal effects 

rs such as excitations to higher order modes at the ends of 
couplers and asymmetric loss in the delay paths, Eq. (4) 
is only an approximate description of the transfer func- 
tion of a real device, and N should be small enough that 
every term in Eq. 4 is meaningful. For our 1.3/1.55jim 

20 WDMs we find that N =3 to 7 is adequate and practical. 
To enhance the stopband rejection > 30 dB, We 
have adopted double- and triple-filtering schemes using 
short chains of N=3 to 7, as illustrated in FIG. 1 (b) and 
(c). FIG. 1 (b) is a double filtered version of (a) and (c) is 

25 a combination of double and triple filtering of four-cou- 
pler and five-coupler chains. For such a multistage filter, 
the transfer function is simply the product of that of each 
stages: 

™ t $> =t $ -t$ • • (8) 

The unwanted light power is thrown away in the extra 
ports instead going to the other output port, and the 
crosstalk is reduced at the expense of rounding the f lat- 
35 ness of the passband. This can be regarded as a further 
optimization of the filter under the constraints of chip 
length and fabrication errors. 

V- System Applications 

40 

In this section, we discuss the potential applications 
of our filters in optical fiber communication systems. 

Our new filter (first, without double or triple filtering) 
has the following properties. First, the two output ports 

45 are complimentary in that the power in the two ports 
always sums up to the input power. Thus the passband 
of one output is the stopband of the other (and for this 
reason we have often used the terms "passband" and 
"stopband" without specifying exactly where they are). 

so Second, the device is symmetric and reciprocal such 
that if the two input ports and simultaneously the two 
output ports are interchanged, or if the input ports are 
interchanged with the output ports, the transmission 
remains the same up to a constant phase factor. Thus 

55 the device only has two distinct transmission states: the 
bar state and the cross state. Third, the device is highly 
directional, i.e., the light propagating in one direction is 
independent of the light propagating in the reverse 
direction. 
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The basic functions of the device are illustrated in 
FIG. 6, where X 1 and are the passband of the bar and 
cross state respectively, or vice versa. (Also, not shown 
here for simplicity, both X 1 and X2 can be groups of 
wavelengths.) As an add-drop filter, it transmits X 1 from 5 
port 1 to port 3, while dropping X 2 into port 4 and add- 
ing the X2 from port 2. Note that the direction of any of 
the signals can be reversed. If X? is absent in pat 1 , the 
device works as a wavelength division multiplexer which 
combines X^ and X 2 in port 3. Similarly, if X 2 is absent /<? 
in port 2, the device works as a wavelength division de- 
multiplexer which separates to port 3 and X2 to port 
4. In the de-multiplexer, if the light in one of the output 
ports is simply discarded, the device works as a simple 
filter. 15 

With double or triple filtering, the situation is more 
complicated. The filter now has more than four ports, 
not all of which are used. However, all the above func- 
tions can still be realized with different configurations 
using the appropriate ports. 20 

Multiplexing and De- multiplexing of the 1.3 and 1.55nm 
Communication Bands 

Most existing fiber communication systems use the 25 
1 .3\im band mainly because fibers have zero dispersion 
and relatively low loss around 1 .3u,m and other compo- 
nents were also first developed for 1 .3gm. Fibers have 
even lower loss away from the dispersion zero around 
1 .55\im, and Er-doped fiber amplifiers (ED FA) are also 30 
readily available in the same wavelength range. There- 
fore, the 1 .55p,m band is expected to coexist with the 
1.3p.m band in many future communication systems, 
and 1 .3/1 .55jim WDM filters will become a key compo- 
nent in such systems. 35 

To accommodate analog-signal transmission and 
future upgrade, 1 .3/1 55fim branching devices used for 
some controlled environments should have passbands 
wider than 1.280-1.335jim and 1 .525-1 .575jim for the 
two output channels respectively, and the crosstalk 40 
should be lower than -50 dB. For other applications, the 
required passbands are even wider (100 nm). Our new 
filters fabricated with the planar waveguide technology 
are ideal candidates for these WDMs, because they 
have the rectangular response required, they are mono- 45 
lithic and reliable, and they can be integrated with other 
components. 

Using the double and triple filtering approach, we 
have designed 1.3/1.55um WDMs meeting the above 
specifications, which correspond to the layouts in FIGs. so 
1(b) and (c). FIG. 1(b) is a double filtered version of (a) 
and FIG. 1 (c) is a combination of double and triple filter- 
ing of four^xjupler and five-coupler chains. The triple fil- 
tering in FIG. 1(c) is only in the 1.3um path to enhance 
the rejection at 1 .55]nm. They have 1 0 or 12 couplers in 55 
series in total and when fabricated with the aforemen- 
tioned doped silica waveguide technology they are 75 
mm long and 0.6 mm wide. For the five-coupler chains 
in FIG. 1(a), (b), and (c), the geometric lengths of the 



couplers are 757, 795, 73, 1209, and 452|im, and the 
geometric path differences are 3.754, 6.256, -5.849, 
and -6.424fim. For the four-coupler chain in FIG. 1(c), 
the geometric lengths of the couplers are 677, 979, 1 99, 
and 1241 jim, and the geometric path differences are 
2.483, 5.733, and -6.055um 

The solid curves in FIG. 7 show our preliminary 
results for the designs shown by the corresponding 
dashed curves. The layout of the three designs is shown 
in FIGs. 1(a), 1(b) and 1(c). For all the designs the 
measured data closely resembles the designed 
response. Wide and flat passbands around both 1.31 
|im and 1 .55 ^im are apparent. Also as designed, the 
stopbands are wide (~100nm) and the transitions are 
sharp. These aspects of our WDMs are better than 
those of any broadband monolithic WDMs previously 
reported. 

Gain Equalization Filters for ED FA Systems 

Er-doped fiber amplifiers (EDFA) have great advan- 
tages over other optical amplifiers used in fiber commu- 
nication systems, but have a highly wavelength- 
dependent gain. For long-range (2= 100km) transmission 
of lightwave signals through optical fibers, EDFAs are 
inserted at spans of every ~50km to compensate the 
attenuation of signals in the fiber. In such a system 
where many EDFAs are cascaded, in order to use the 
full bandwidth of the EDFAs (1.53-1.56 ^im), a gain 
equalization filter must be used along with each EDFA 
to flatten the overall system gain. The response of these 
filters is roughly the inverse of the gain of the EDFA, and 
has a peak at 1 .538 fim with asymmetric wings. Despite 
the irregular shape required, these filters can be readily 
designed with our Fourier expansion approach and fab- 
ricated using the planar waveguide technology. 

The solid curve in FIG. 8 shows a designed EDFA 
gain equalization f iher using a seven-coupler chain. The 
circles represent the required filter response, optimized 
for maximum end-to-end flatness over a 30 nm band- 
width, for a system with 40 km spans, 2000 km total 
length, using two-stage EDFA (type E200S) with dual 
90 mW 980 nm pumps. Because of the sharp features 
around 1.54 ^im, a 8-ratio of approximately 1/-2/-2/-2/- 
4/-4 is found to be better than 1/2/2/2/2/2. The geomet- 
ric lengths of couplers are 1002, 861, 571, 1112, 1367, 
842, 1180fim, and the geometric path differences are 
13.56, -27.09, -26.93, -26.80, -53.16, -53.70^m. The 
design is also relatively stable against overall errors in 
the coupling length. 

Claims 

1. A monolithic optical waveguide filter comprising a 
pair of optical waveguides forming a sequence of 
N>2 couplers alternating with a plurality of N-1 
delay paths between said couplers, wherein the 
effective length of one of said couplers differs by 
more than 5% from the effective length of another 
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coupler and the differential delay provided by one of 
said delay paths differs by more than 10% from the 
differential delay provided by another of said delay 
paths. 

5 

2. A monolithic optical waveguide filter according to 
claim 1 wherein the differential delay provided by at 
least one of said delay paths is opposite in sign 
from the differential delay provided by another of 
said delay paths. 10 

3. A monolithic optical waveguide filter comprising a 
pair of optical waveguides forming a sequence of 
N>2 couplers alternating with a plurality of N-1 dif- 
ferential delay paths between said couplers, is 
thereby providing a plurality of optical paths 
between an input at one end of said sequence and 

an output at the other end, wherein each path of 
said plurality presents at said output an optical sig- 
nal whose amplitude and phase form a term of a 20 
Fourier series comprising the filter transmission 
function. 

4. A monolithic optical waveguide filter according to 
claim 3 wherein each path presents a term in an 25 
odd Fourier series. 

5. A monolithic optical waveguide filter comprising a 
pair of optical waveguides forming a sequence of 
N>2 couplers alternating with a plurality of N-1 30 
delay paths between said couplers, wherein said 
delay paths provide, in any order, within ± 25%, one 
differential delay of proportion ± 1 and one or more 
differential delays of proportion ± 2 or ± 4 or ± 8. 

35 

6. A monolithic optical waveguide filter according to 
claim 5 wherein the effective length of one of said 
couplers differs by more than 5% from the effective 
length of another coupler. 

40 

7. A monolithic optical waveguide filter according to 
claim 5 wherein the differential delay provided by 
one of said delay paths is opposite in sign from the 
differential delay provided by another of said delay 
paths. 46 

8. A monolithic optical waveguide filter comprising a 
pair of optical waveguides forming a sequence of 
N>2 couplers alternating with a plurality of N-1 
delay paths between said couplers, wherein said so 
delay paths provide, in any order, within ± 25%, two 
differential delays of proportion ± 1 and one or more 
differential delays of proportion ± 2, or ± 4 or ± 8. 

9. A monolithic optical waveguide filter according to ss 
claim 8 wherein the effective length of one of said 
couplers differs by more than 5% from the effective 
length of another coupler. 



10. A monolithic optical waveguide filter according to 
claim 8 wherein the differential delay provided by 
one of said delay paths is opposite in sign from the 
differential delay provided by another of said delay 
paths. 

1 1 . A filter according to claim 1 or 2 or 3 or 4 or 5 or 8 
wherein 3 5 N 5 7. 

12. A monolithic optical waveguide filter according to 
claim 1 further comprising one or more additional 
optical filters optically coupled to said pair of optical 
waveguides to form a filter network. 

13. A monolithic optical waveguide filter according to 
claim 4 or 1 2 wherein said light from said paths pro- 
vides a rectangular filter response. 

14. A monolithic optical waveguide filter according to 
claim 4 wherein said light from said paths provides 
a rectangular filter response for separating 1.3 
micrometer wavelength light from 1.55 micrometer 
wavelength light. 

15. A monolithic optical waveguide filter according to 
claim 4 wherein said light from said paths provides 
a filter response inverse to the gain spectrum of an 
erbium-doped fiber amplifier. 
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